Stochastic modeling of regulation of gene expression by multiple small RNAs 
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A wealth of new research has highlighted the critical roles of small RNAs (sRNAs) in diverse 
processes such as quorum sensing and cellular responses to stress. The pathways controlling these 
processes often have a central motif comprising of a master regulator protein whose expression is 
controlled by multiple sRNAs. However, the regulation of stochastic gene expression of a single 
target gene by multiple sRNAs is currently not well understood. To address this issue, we analyze 
a stochastic model of regulation of gene expression by multiple sRNAs. For this model, we derive 
exact analytic results for the regulated protein distribution including compact expressions for its 
mean and variance. The derived results provide novel insights into the roles of multiple sRNAs in 
fine-tuning the noise in gene expression. In particular, we show that, in contrast to regulation by 
a single sRNA, multiple sRNAs provide a mechanism for independently controlling the mean and 
variance of the regulated protein distribution. 
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Small non-coding RNAs (sRNAs) are known to play a 
central role in diverse cellular pathways that bring about 
global changes in gene expression [1]. In several cases, 
such global changes are coordinated by a master regu- 
lator protein whose expression is controlled by multiple 
sRN As [2, ^] . Examples include regulation of the master 
regulator in bacterial quorum-sensing pathways by mul- 
tiple sRNAs [i^ and regulation of the alternative sigma 
factor by four distinct sRNAs, each of which responds 
to different environmental stresses [2] . Despite its impor- 
tance in coordinating cellular stress responses and related 
processes, the role of multiple sRNAs in regulating the 
expression of a single target gene is not yet well under- 
stood [3]. In this work, we address this issue by analyzing 
a stochastic model that elucidates potential functional 
roles for this widely-occurring regulatory motif. 

Regulation of gene expression by sRNAs is a post- 
transcriptional process: sRNAs can bind to messenger 
RNAs (mRNAs) and control protein production by al- 
tering mRNA stability or by regulating translational ef- 
ficiency [l|. The intrinsic stochasticity of the underlying 
biochemical reactions can produce significant variations 
('noise') in gene expression among individual cells in iso- 
genic populations [5|49|. Although noise in gene expres- 
sion can have deleterious effects in some cases and thus 
needs to be limited; in other cases such noise is utilized 
and indeed required by the cell e.g . for processes leading 
to cell-fate determination [iq|, 11 |. Furthermore, it has 
been argued that noise in gene expression could be advan- 
tageous under conditions of high stress, since variability 
in a population provides a bet-hedging strategy that can 
enable survival [ll, E^. Regulation of the noise in gene 
expression is thus essential for the proper functioning of 
several cellular processes. Since sRNAs regulate critical 
cellular processes, understanding their role in fine-tuning 
the noise in gene expression is of fundamental importance 

lil. 



A quantitative understanding of the cellular functions 
of sRNAs is aided by the development of models, which 
can often produce insights that guide future experiments. 
In recent research, several models for regulation by sR- 
NAs have been developed. Since many sRNAs are known 
to repress gene expression, most previous models have fo- 
cused on re gul ation by irreversible stoichiometric degra- 
dation 14l4l7|. However, sRNAs can affect not only 



mRNA degradation rates but also protein production 
rates and the corresponding biochemical reactions are, 
in general, reversible [l], [l8|. Furthermore, not all sR- 
NAs repress gene expression; there are sRNAs which are 
known to activate gene expression and even some which 
can switch from activating to repressing in response to 
cellular signals pj|. To quantify the corresponding 
effects on stochastic gene expression, a general model 
which includes the different mechanisms of sRNA-based 
regulation needs to be analyzed. Such a model, for the 
case of a single sRNA regulator, has been developed in re- 
cent work [20] . Analysis of this model and its extension to 
multiple sRNAs thus provides a means of addressing out- 
standing questions about the impact of different modes 
of regulation by sRNAs on the noise in gene expression. 

In this work, we generalize our previous model [20] to 
analyze the case of multiple sRNAs regulating a single 
mRNA target. Specifically, we derive exact analytic ex- 
pressions for the generating function of a protein burst 
distribution resulting from the regulation of a single tar- 
get by an arbitrary number of sRNAs. Using this expres- 
sion, we obtain compact analytic expressions for both the 
mean and variance of the regulated protein distribution. 
We first analyze these results for the case of a single reg- 
ulator and derive insights into different modes of regula- 
tion by sRNAs. These results are then contrasted with 
features unique to regulation by multiple sRNAs. In par- 
ticular, we show that, in contrast to regulation by a single 
sRNA, regulation by multiple sRNAs provides the cell 
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FIG. 1. Schematic illustration of regulation of gene expression 
by multiple sRNAs. In the full reaction scheme, there are N 
different regulators and the kinetic scheme is shown for the 
i^^ sRNA regulator. The association and dissociation rates for 
binding to the mRNA are denoted by ai and /Si respectively. 
Association results in a complex which produces proteins with 
rate kp. and is degraded with rate /j.^. Note that for the 
mRNA to transition from one complex to another, it must first 
return to its unbound state before forming a new complex. 



with a mechanism to independently fine-tune both the 
mean and variance of the regulated protein distribution. 

We begin by considering protein production from a 
single mRNA regulated by N independent sRNAs. The 
corresponding reaction scheme is shown in Figure 1. The 
mRNA has TV + 1 possible states, with the states i = 
1, .., A/" denoting mRNA bound to the i^^ sRNA regulator 
to form complex i. For notational simplicity, we denote 
the unbound mRNA state as complex 0. An unbound 
mRNA forms complex i with rate ai and the complex can 
either dissociate with a rate decay with a rate fia^ or 
initiate protein production with a rate kp.. We assume 
the sRNA regulators are present in large amounts such 
that fluctuations in their concentration can be ignored; 
correspondingly the rates are taken to be constant. 

The distribution of proteins produced from a single 
mRNA (interacting with N sRNAs) before it decays is 
denoted as the protein burst distribution P5^Ar(n). We 
further define the function fi{n,t) which denotes the 
probability that n proteins have been produced and the 
mRNA is in state i at time t. Correspondingly, the pro- 
tein burst distribution, P5^Ar(n), can be obtained from 
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Pb,N{n) = ^ fi(n,t)iicidt. 



(1) 



The time-evolution of the probabilities fi{n^t) is gov- 
erned by the Master equation 
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The initial condition corresponds to creation of a sin- 
gle unbound mRNA and no proteins in the system at 
time t = 0, i.e. /o(0,0) = 1. The procedure outlined 
in [20] can now be applied to obtain the burst distri- 
bution or more precisely the corresponding generating 
function Gb,N{z) = ^^-^b,N{n). Specifically, defining 
Fi{z^ s) = e~^*fi{n^ t) dt^ we obtain the gener- 

ating function using 
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Gb{z) = \im'Y(^iJ.c,Fi{z, s)^ 



(3) 
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To present the results in a compact form, it is convenient 
to define the dimensionless variables = ^ and 



(^) for i > and = ^ for z > 0. Now, 

by setting ujq = 1 and = we further define the 'weight 
functions' uJi{z) = Ui . Note that is 

the generating function of a geometric distribution with 
mean ^i. 

Using the above definitions, we obtain the following ex- 
act expression for the generating function of the protein 
burst distribution 
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For A'' = 0, i.e. the unregulated case, the generating func- 
tion reduces to 
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l + no(l-z) 



(5) 



in agreement with previous work showing that the pro- 
tein burst distribution is a geometric distribution with 
mean no = jf^ 0, HH . Eq. (j4]) provides the general- 
ization of this result for the case of N sRNA regulators 
using the weight functions 0Ji{z). 

An important mechanism of regulation by sRNAs cor- 
responds to the case that sRNA binding prevents ribo- 
some access and thus blocks translation. For the case 
that all the regulators act to fully repress translation, 
i.e. kp. = for z > 0, we have uji{z) = uji. Correspond- 
ingly, the generating function reduces to the formula for 
the unregulated case (Eq. (|5j)) with a renormalized mean 
given by ^J}^ — . This interesting observation indicates 

that regulation by sRNAs which function by fully repress- 
ing translation is reversible: for arbitrary concentrations 
of the sRNA regulator, by appropriately adjusting the 
parameter kp^ (e.g. by adjusting the concentration of ri- 
bosomes), the regulated protein distribution in the pres- 
ence of sRNAs can be made identical to the distribution 
for the unregulated case (prior to introduction of the sR- 
NAs). 

For the general case, the generating function can be 
recast in a form that shows that the protein burst distri- 
bution is a mixture of A/" + 1 geometric distributions [20i] . 



However, the corresponding expression, even for the case 
of = 2, is too complex to be reproduced here. On the 
other hand, using Eq. (H]), compact analytic expressions 
for the mean, 77,5^, and squared coefficient of variance, 
o'^n/'^In derived. The mean (scaled by the un- 

regulated mean) is given by 
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(6) 



and the noise strength (squared coefficient of variance) is 
given by 
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where 



Y,i=o {rij - no) 



Qn 



Note that the signs of Fn and Qn characterize the im- 
pact of the sRNAs on the regulated protein distribu- 
tion. Specifically, the unregulated case has mean n^^; 
thus Fn < corresponds to repression whereas Fn > 
corresponds to activation. Similarly, an unregulated pro- 
tein burst distribution with mean n^^ has a squared co- 
efficient of variance 1 + l/n^^; thus, when Qat < we 
have noise reduction whereas Q at > corresponds to in- 
creased noise strength (relative to an unregulated burst 
distribution with the same mean). 

We now focus on using Eq. (|6j) and Eq. ([7]), to eluci- 
date interesting features for the case of regulation by a 
single sRNA, i.e. N = 1. Note that all of the variables 
in the expressions for the mean and noise strength are 
always positive (or zero) except for the term (ni — no). 
Thus, the sign of Fi and Qi is determined completely 
by Aio = ni — no- When Aio > both the mean and 
the noise strength are higher than their unregulated val- 
ues. Similarly, when Aio < both the mean and the 
noise strength are lower than the corresponding unreg- 
ulated values (except for the case = for which the 
noise strength is identical to an unregulated distribution 
with the same mean). In either case, we note that, for 
a single sRNA regulator, there is a coupling between the 
mean and noise strength of the regulated burst distri- 
bution such that both cannot be tuned independently, 
e.g. a decrease in the mean cannot be associated with an 
increase in the noise strength. 

In contrast to the case of regulation by a single sRNA, 
in the case of regulation by multiple sRNAs, the mean 
and noise of the protein distribution can be tuned in- 
dependently. The deviation of the mean from its un- 
regulated value depends solely upon terms of the form 
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FIG. 2. Contour plots for the percent change in the mean 
and noise strength of a two regulator system from the corre- 
sponding unregulated values as a function of /cp^ and kp^ . (a) 



Mean: Plot of f{kp-^, kp 



100%. Note that along 



the contour f{kp^^kp^) = —20% the noise strength changes 
from less than -5% to over 70%. (b) Noise Strength: Plot 
of g(kp-^^kp^) — Y^fni, ' Note that g(kp-^^kp^) con- 

tains contours that sweep out a large portion of the plotted 
(A^Pi, kp^) state space. By proportionally changing the kp val- 
ues corresponding to the two regulators, the noise strength 
can be varied while maintaining the same mean value. The 
parameters used were kp^ — 50, /ico = I5 Mci = 4.5, /ic2 — 4.5, 
/3i = 1, /32 = 0.5, ai = 2 and 0^2 2. 



AiO = — uq. On the other hand, considering the gen- 
eral form of the noise strength for N > 1, we have terms 
of the form Aij = (^^ — ij){ni — nj) that contribute to 
the deviation from the corresponding unregulated value. 
Thus, for appropriately chosen parameters, two sRNAs 
can be used to tune both the mean and variance of the 
regulated protein distribution as discussed below. 

Consider the case of regulation by 2 sRNAs that are 
maintained at some fixed cellular concentrations. A new 
mRNA target for these sRNAs can arise from the evo- 
lution of appropriate sRNA binding sites on the mRNA 
sequence. For the new target, we assume that the param- 
eters kp^ and can be tuned based on changes in the 
sequence and location of the sRNA binding sites. The 
corresponding variation in the mean and noise strength 
is shown in Fig. 2. Note that by maintaining a linear re- 
lationship between kp^ and kp^ , the mean of the regulated 
protein distribution can be left unchanged; however, the 
noise strength can be tuned over a large range. For ex- 
ample, for some choices of the parameters, the mean can 
be fixed and the noise strength can be varied by over 
100% relative to the unregulated distribution (see Fig. 
2). In this context, it is interesting to note that it has 
been observed that several sRNAs have a minimal effect 
on the mean levels of their regulatory targets. For such 
targets, sRNAs could be functioning primarily as modu- 
lators of noise while giving rise to only a minimal change 
in mean levels due to regulation [11|. Our results pro- 
vide quantitative insight into how such regulation can be 
implemented using multiple sRNA regulators. 
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FIG. 3. Contour plot of the percent change in noise strength 
of a two regulator pathway from its corresponding unreg- 
ulated value as a function of ai and a2, i.e. f{ai,a2) = 
— • 100%. The change in mean from the unregulated 

to regulated pathways, Fn, is positive below and negative 
above the line ai = a2- The parameters used were kpQ — 50, 
kp^ = 200, kp^ = 72.5, /Xco = 1, Mci = 2.725, /Xci = 2.725, 
/3i 0.15 and /32 0.15. 

The results obtained also illustrate how changing 
sRNA concentrations can be used to modulate the noise 
in gene expression. For our model, changes in the concen- 
tration of the sRNA regulators effectively alter the bind- 
ing rates {ai) to the mRNA. From Eq. (|6j), we see that for 
2 regulators, by choosing one of the regulators to be a re- 
pressor and the other to be an activator, the mean of the 
regulated protein distribution can be increased {Fn > 0) 
or decreased {Fn < 0) by adjusting the relative concen- 
trations of the two regulators. Furthermore, by changing 
the concentrations of the regulators such that their rel- 
ative concentration is fixed, the mean of the regulated 
protein distribution is left unchanged, whereas the vari- 
ance can be tuned over a range of values. This insight 
is particularly relevant, given that noise can be advanta- 
geous to a cell, especially in response to stress. 

Finally, we note that while the above analysis focuses 
on the burst distribution from a single mRNA, the results 
are readily generalized to the case of arbitrary mRNA 
burst distributions, assuming independent contributions 
from each mRNA [22]. Assuming independent bursts and 
using equations derived in recent work [22j], we can use 
the derived results to obtain corresponding expressions 
for the mean and variance of the steady-state protein 
distribution. A detailed analysis of the connection to 
steady- state distributions will be presented elsewhere. 

Modulation of gene expression in pathways containing 
multiple sRNA regulators could provide a multitude of 
selective advantages in specific contexts and our model 
provides a means of gaining insight into these processes. 
In particular, our analysis shows how by adjusting the 



concentrations of multiple sRNA regulators, a cell can 
initiate finely tuned responses to external stimuli. This 
could explain the ubiquity of sRNA regulators in cellular 
stress response, a process for which gene expression 
noise is known to be critical. In a broader context, the 
fine control afforded to cells via regulation by multiple 
sRNAs could be useful in signal integration when mul- 
tiple environmental stimuli are present simultaneously. 
This work provides a basis for future work investigating 
such signal integration and response in more complex 
pathways; thus opening new avenues for understanding 
and modeling stochastic gene expression in a wider class 
of regulatory networks. 
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